Hybrid stars within a covariant, nonlocal chiral quark model 
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We present a hybrid equation of state (EoS) for dense matter in which a nuclear matter phase is 
described within the Dirac-Brueckner-Hartree-Fock (DBHF) approach and a two-flavor quark matter 
phase is modelled according to a recently developed covariant, nonlocal chiral quark model. We show 
that modern observational constraints for compact star masses (M ~ 2Mq) can be satisfied when a 
small vector- like four quark interaction is taken into account. The corresponding isospin symmetric 
EoS is consistent with flow data analyses of heavy ion collisions and points to a deconfmement 
transition at about 0.55 fm -3 . 

PACS numbers: 04.40.Dg, 12.38.Mh, 26.60.+C, 97.60.Jd 



I. INTRODUCTION 



Understanding the properties of matter at moderate and high densities is required, e.g., to explain the astrophysical 
phenomena which accompany the birth of compact stars in supernova explosions, and the further evolution processes 
of cooling, spin-down, accretion, merging with companion stars, etc., which lead to effects accessible to observation. 

Nowadays, one of the questions in the focus of discussions is the possibility of a phase transition to deconfined quark 
matter in the stellar cores. There is no doubt that deconfmement of quarks shall occur at sufficiently high densities, in 
accordance with the asymptotic freedom of QCD [1]. However, it is not clear a priori whether the critical density for 
deconfmement is low enough to be reached in the cores of neutron stars [2|. Even if this condition can be fulfilled, the 
properties of dense quark matter in the vicinity of the deconfmement transition might be too similar to those of dense 
hadronic matter to result in clearly distinguishable signals. On the other hand, there have been a number of interesting 
suggestions about how the occurrence of quark matter in neutron stars could manifest itself and possibly contribute 
to the resolution of puzzling observations. For the sake of illustration we remind on the timing behavior of pulsar 
spin-down frequency clustering [i[ or population clustering [1, @ of accreting compact stars, which are based on 
a softening of the equation of state and therefore a compactification as well as a reduction of the maximum allowable 
stellar masses. Since deconfined quark matter is rather stiff when compared to hyperonic matter, the problem with 
rather low maximum masses of compact stars with hyperonic interior could be solved by the occurrence of a quark 
matter core ■ Similarly, the description of the compact star cooling evolution with superconducting quark matter 
interior seems favorable over a purely hadronic modelling [8j. A deconfmement transition during the protoneutron 
star (PNS) evolution offers a mechanism to explain the gamma-ray burst energy release of the order of 100 bethe 
(=10 53 erg) H, 0, Hi] . The nucleation timescales for a quark matter phase transition could explain the time delay 
statistics of GRB subpulse structure [H, EH- In the presence of a strong magnetic field, neutrino propagation in hot, 
superconducting quark matter can become collimated (beaming) and asymmetric, thus explaining a resulting kick 
velocity for the PNS [3. 
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The recent progress in compact star observations justifies a reinvestigation of the issue of hybrid stars with quark 
matter cores and theoretical aspects of dense matter properties. In particular, the high mass of M = 2.1 ± 0.2 Mq 
for the pulsar J0751+1807 in a neutron star-white dwarf binary system [l5[ and the large radius of R > 12 km 
for the isolated neutron star RX J1856. 5-3754 (shorthand: RX J1856) [16] point to a stiff equation of state at high 
densities. Measurements of high masses are also reported for c omp act stars in low-mass X-ray binaries (LMXBs) as, 
e.g., M — 2.0 ±0.1 Mq for the compact object in 4U 1636-536 [17| . With data of this kind, new st ring ent constraints 
on the equation of state of st rong ly interacting matter at high densities have been formulated (see [181 ] and references 
therein) . It has been argued [161 Il9j that deconfined quark matter cannot exist in the centers of compact stars with 
masses and radii as reported for these objects. In view of recent works on the quark matter EoS, however, this claim 
app ears to be premature [2(J. It has been demonstrated within the Nambu-Jona-Lasinio model for quark matter 
[2l| that the inclusion of a diquark condensate (leading to color superconductivity) together with a vector meson 
condensate, not only provides a more elaborate description of the EoS but also allows to describe the phenomenology 
of hybrid stars in excellent accordance with the above mentioned new mass and mass-radius constraints. The two 
mechanisms at work are: (1) the lowering of the phase transition density due to the diquark condensate, so that 
already typical neutron stars with masses in the range 1.1 — 1.5 Mq can be hybrid stars with extended quark matter 
cores; (2) the stiffening of the EoS due to the vector mean field, which implies an increase of the maximum accessible 
masses of star configurations up to ~ 2 Mq. 

Another lesson to be learnt from NJL model studies is that at low temperatures there is a sequential deconfinement: 
strange quarks occur only at densities well above the deconfinement of light quarks [22], [23J, |24j, |25|, |26j, [27J. If 
those densities could be reached in a compact star, the corresponding strange quark matter cores would be in a 
superconducting CFL phase, which renders the hybrid star configuration mechanically unstable [2l], [28|, [2j|. In the 
present work we will restrict ourselves to the discussion of the two-flavor case, applying the more elaborate formalism 
of a recently developed nonlocal, covariant chiral quark model (30l. I3lll32l] . Moreover, we consider a generalized version 
of this model, including an isoscalar vector meson current which, in the same way as in the case of the Walecka model 
for nuclear matter, leads to a stiffening of the quark matter EoS. 

Although we understand hadrons as bound states of quarks, there does not yet exist a unified approach which 
accurately describes the thermodynamics of the transition from nuclear to quark matter. Therefore, we apply a 
two-phase description by performing a Maxwell construction describing the transition from a nuclear matter EoS to 
the quark matter EoS. The nuclear equation of state to be considered here results from calculations within the DBHF 
approach [33j and has already been applied for the description of compact stars before [la , \2l\ . Since this EoS is 
rather stiff, several modern compact star observations have been well reproduced. In particular, within this approach 
one is able to obtain large neutron star masses (M max — 2.33 Mq) and radii (R = 12 — 13 km for typical neutron 
stars). On the other hand, the DBHF description is not so well suited to reproduce the results obtained from elliptic 
flow data in symmetric nuclear matter (SNM): from these results, it can be seen that the DBHF EoS tends to be too 
stiff beyond densities of 3 times the saturation density n sat = 0.16 fm -3 jl8j. We show here that this problem can be 
solved by a phase transition to quark matter while simultaneously fulfilling the constraints on the behavior of dense 
matter for the case of hybrid stars. 

The article is organized as follows: a brief description of the DBHF approach and the quark matter model used 
here is given in Sects. II. A and II. B respectively. In Sect. Ill we show our numerical results, comparing them with 
present empirical constraints on the behavior of dense matter under constraints for conditions in neutron stars and 
heavy ion collisions. Our conclusions are stated in Sect. IV. Finally, in the Appendix we provide some details of the 
quark matter model. 

II. THEORETICAL FORMALISM 

As mentioned in the Introduction, in the present calculation we use a two phase description to account for the 
transition from a nuclear matter EoS to a quark matter EoS. In the following two subsections we briefly discuss the 
theoretical approaches considered here to describe each of these phases. 

A. Nuclear matter equation of state 

For the description of the nuclear matter equation of state we rely on the relativistic Dirac-Brueckner-Hartree-Fock 
(DBHF) approach where the nucleon inside the medium is dressed by the self-energy X. This self-energy is obtained 
from the Bethe-Salpeter equation for the nucleon-nucleon T-matrix in the ladder approximation, using the Bonn-A 
potential in the interaction kernel (34[. 
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We employ a parameterization of the resulting EoS based on a parabolic dependence of the energy per nucleon on 
the asymmetry parameter a = 1 — 2x, given in the form 

E(n,a) = E (n) + a 2 E s (n) , (1) 

where x = n p /n is the proton fraction, Eq(ti) is the energy per nucleon in SNM, and Es{n) is the (a)symmetry energy. 
Both contributions Eo(n) and Es(n) have been extracted from DBHF calculations for the cases a = and a = 1, 
respectively. The parabolic interpolation in Eq. (fT]) has been widely used in the literature, see e.g. Refs. flRl35l|. and 
proves to be an excellent parameterization of the asymmetry dependence for the purpose of the present study. An 
exact reproduction of a given EoS might require higher order terms, which have been neglected here. The advantage 
of the parabolic interpolation lies in the fact that all zero temperature equations of state for neutron star matter 
(NSM) can be derived by applying simple thermodynamic identities [36j . In particular, we obtain 

es(n,a) = nE(n,a) , (2) 
d 

P B {n,a) = n 2 —E(n,a) , (3) 

Hn,p{n,a) = (^l + n^^ E (n) ~ (^a 2 T^a~ a 2 n^j E s (n) 7 (4) 

for the baryonic energy density e(n), the pressure P{n), and the chemical potentials of neutrons /x„ (upper sign) and 
protons Hp (lower sign), respectively. 

B. Quark matter equation of state 

Early reviews on the treatment of quark matter within the NJL model as a chiral quark model with a local current- 
current-type interaction can be found, e.g., in Refs. [37l. l38l. l39j . In order to include color superconductivity, techniques 
were needed which are described, e.g., in Refs. j40Ll41j. 

We use here a generalization of these approaches to a nonlocal chiral quark model which includes scalar and vector 
quark- antiquark interactions and anti-triplet scalar diquark interactions. The corresponding effective Euclidean action 
in the case of two light flavors is given by 

l x {^(x) ( -^ + m J^)-^(4^ • (5) 

Here m c is the current quark mass, which is assumed to be equal for u and d quarks, whereas the currents js.D,v{x) 
are given by nonlocal operators based on a separable approximation to the effective one gluon exchange model (OGE) 
of QCD. These currents read 



f s (x) = J d^z g(z)Tjj(x + -) T f 1>(x--), 
fb( x ) = J d 4 z g(z) 4>c(x + ~) «75T- 2 A a ip(x - |) , 

3v( x ) = J d * z 9( z ) + |) 7m V>(z - |) , (6) 

where we have defined ipc(%) = 7274 ^ T {x) and Tf = (ljij^f), while f and A a , with a = 2,5,7, stand for Pauli 
and Gell-Mann matrices acting on flavor and color spaces, respectively (notice that 7 M = (7,74) are Euclidean Dirac 
matrices). The functions g(z) in Eqs. ([6]) are nonlocal covariant form factors characterizing the effective quark 
interaction [32I ]. 

The effective action in Eq. |5J might arise via Fierz rearrangement from some underlying more fundamental in- 
teractions, and is understood to be used — at the mean field level — in the Hartree approximation. In general, the 
ratios of coupling constants H/Gs, Gy/Gs would be determined by this microscopic couplings; for example, OGE 
interactions in the vacuum lead to H/Gs — 0.75 and Gy/Gs = 0.5. However, since the precise derivation of effective 
couplings from QCD is not known, there is a large theoretical uncertainty in these ratios. Details of the values used 
in the present work will be given below. 

We proceed by considering a bosonized version of this quark model, in which scalar, vector and diquark fields 
are introduced. Moreover, we expand these fields around their respective mean field values, keeping the lowest order 
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contribution to the thermodynamic quantities. The only nonvanishing mean field values in the scalar and vector sectors 
correspond to isospin zero fields, a and Co respectively, while in the diquark sector, owing to the color symmetry, one 
can rotate in color space to fix A5 = A 7 = 0, A 2 = A. 

Now we consider the Euclidean action at zero temperature and finite baryon chemical potential \ib ■ Introducing dif- 
ferent chemical potentials /i/ c for each flavor and color, the corresponding mean field grand canonical thermodynamic 
potential per unit volume can be written as 

where the inverse propagator S^ 1 is a 48 x 48 matrix in Dirac, flavor, color and Nambu-Gorkov spaces [13, [H[. Its 
explicit regularized form ^^eg) ^ s gi ven m the Appendix, together with further details about the model. The mean 
field values a, A and Co are obtained from the coupled equations 

r]Q MFA r!Q MFA HQ MFA 

^A- = < = < ^F- = ■ ( g ) 

c(A da aui 

In principle one has six different quark chemical potentials, corresponding to quark flavors u and d and quark 
colors r, g and b. However, there is a residual color symmetry (say, between red and green colors) arising from the 
direction of A in color space. Moreover, if we require the system to be in chemical equilibrium, it can be seen that 
chemical potentials are not independent from each other. In general, it is shown that all fif c can be written in terms 
of three independent quantities: the baryonic chemical potential fj,g, a quark electric chemical potential fiq and a 
color chemical potential /^s- The corresponding relations read 

HB 2 1 ji B| 2 2 

Mur — Pug — — + gMQ, + gM8 ! M«6 — ~ + gMQg _ gM8 , 

MB 1 1 MB 1 2 
Mdr - Mdg = ~^ gMQ, + gM8 , Mdh - ~g gMQ, - gM8 • (yj 

The chemical potential MQ,; which distinguishes between up and down quarks, as well as the color chemical potential 
M8j which has to be introduced to ensure color neutrality [42L l43l. l44l|. vanish for an isospin symmetric quark matter 
system. Thus, in this case, the corresponding EoS can be obtained after calculating the mean field values ct, Q and A 
from Eqs. ©. 

Now, if we want to describe the behavior of quark matter in the core of neutron stars, in addition to quark matter 
we have to take into account the presence of electrons and muons. Thus, treating leptons as a free relativistic Fermi 
gas, the total pressure of the quark matter + lepton system is given by 

p = -n^j-n', (10) 

where Vl l is the thermodynamical potential per unit volume for a gas of noninteracting electrons and muons (see 
Appendix A). In addition, it is necessary to take into account that quark matter has to be in /3 equilibrium with 
electrons and muons through the (3 decay reactions 

d~ >u + l + Di, u + l^d+vt, (11) 

for / = e,/U. Thus, assuming that (anti)neutrinos escape from the stellar core, we have an additional relation between 
fermion chemical potentials, namely 

[Uc - Hue = — MQ 9 = Mi ( 12 ) 

for c = r, 5, b, fi e = /i M = 

Finally, in the core of neutron stars we also require the system to be electric and color charge neutral, hence the 
number of independent chemical potentials reduces further. Indeed, fii and fig, get fixed by the condition that charge 
and color densities vanish, 

PQtot = PQ q ~ ^2 Pi = X] ( \ Puc ~ \ Pdc ) ~ P l = ' 



= ~m ( p f r + p f9 ~ 2 Pf b "> = > ( 13 ) 

f=u,d 
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where the expressions for the lepton densities pi and the quark densities pf c can be found in the Appendix. Note 
that the set of color chemical potentials required to ensure color neutrality of the system depends on the choice of the 
orientation of the diquark condensate orientation in color space [45] . For the standard choice employed in the present 
work, /is is sufficient. In summary, in the case of neutron star quark matter, for each value of ps one can find the 
values of A, a, u>, fn and ps by solving Eqs. ([8]), supplemented by Eqs. (fl"2|) and (fl"3|) . This allows to obtain the quark 
matter EoS in the thermodynamic region we are interested in. 

III. NUMERICAL RESULTS AND DISCUSSION 

In this section we present our numerical results, showing the behavior of both isospin symmetric hadronic matter and 
neutral hadronic matter for finite baryochemical potential. The nuclear-to-quark matter phase transition is treated 
in the traditional way, following a two-phase scheme in which the nuclear and quark matter phases are described by 
the theoretical approaches presented in the previous section. 

As stated, the nuclear matter phase is described according to the DBHF approach, in which the nucleon self-energy 
is calculated from the Bethe-Salpeter equation considering a Bonn-A potential in the interaction kernel. Regardin g th e 
quark matter model, we note first that, according to previous analyses carried out within nonlocal scenarios [30Ll3ll| . 
the results are not expected to show a strong qualitative dependence on the shape of the nonlocal form factors. Thus 
we will consider (in momentum space) a simple and well-behaved Gaussian function, 

g(p 2 ) = exp(-p 2 /A 2 ) , (14) 

where A is a free parameter of the model, playing the role of an ultraviolet cut-off. The value of A, as well as the 
values of the free model parameters Gs and m c , can be fixed from low energy phenomenology. Here we have chosen 
these input parameters so as to reproduce the empirical values for the pion mass m w = 139 MeV and decay constant 
fir = 92.4 MeV, and to obtain a phcnomcnologically reasonable value for the chiral condensate at vanishing ps, namely 
(0\qq\0)^ 3 = -230 MeV. In this way we obtain m c = 6.49 MeV, G s = 0.515 x 10" 4 MeV" 2 and A = 678 MeV [H. 

The values of the scalar diquark coupling and the isoscalar vector coupling H and Gy, or equivalently the di- 
mensionless ratios h = H/G,s and g = Gv/Gs, are considered here as parameters to be chosen in accordance with 
phcnomenological constraints from flow data analyses of heavy-ion collisions and the new mass and mass-radius con- 
straints from compact star observations. In accordance with our previous investigation of the phase diagram of neutral 
quark matter in Ref. [321 ] we choose two typical values h — 0.70 and h = 0.74 for the diquark/scalar coupling ratio, 
and analyze the phase transition features for different values of g so as to obtain an acceptable transition density. 

In Fig. [1] we show the curves for the pressure as function of the baryochemical potential, for both nuclear and quark 
matter phases, in the case of isospin symmetric matter. Left and right panels correspond to h — 0.70 and h = 0.74 
respectively, and in each case three values of the scaled vector coupling g have been chosen. As can be seen from 
the figures, within these ranges of g both nuclear and quark matter EoS behave similarly in the relevant domain of 
baryochemical potentials, therefore the value of the critical \ib strongly depends on the value of the vector coupling. 
In any case, for both h = 0.70 and h = 0.74 it is possible to tune the value of g so that the softening due to the 
deconfinement transition occurs at a density of about 0.55 fm -3 . At this point the DBHF EoS becomes too stiff, 
being unable to fulfill the flow constraint [47| . This situation is sketched in Fig. [21 where we show the area allowed by 
the flow constraint and the corresponding DBHF curve in the density-pressure plane. For h = 0.70 and h = 0.74, the 
values g = 0.05 and g = 0.07 respectively (dotted lines in the P — fin curves of Fig.[T]) lead to nuclear-to-quark matter 
phase transitions such that the stiffness of the pressure curve is softened and the flow constraint can be satisfied. On 
the other hand, values of g outside the range considered in Fig. [T] would lead to either a too early phase transition or 
to a situation in which there is no transition to quark matter at all. We emphasize that both the nuclear and quark 
matter EsoS behave almost identically over a wide range of the baryochemical potential in Fig. [TJ This results in an 
almost direct crossover transition with a small or even negligible density jump at the transition density, as shown in 
Fig. [2 As a further consequence, the transition density strongly depends on small changes of g. 

Finally, let us take into account the hadronic matter models leading to the curves in Fig. [2l considering now 
the description of neutral hadronic matter in the interior of compact stars. In Fig. [31 we show the sequences of 
compact star configurations obtained as solutions of the Tolman-Oppenheimer-Volkoff equations of general-relativistic 
hydrodynamic stability, for selfgravitating dense hadronic matter described by these theoretical models. The results 
quoted in Fig.[3jrepresent the main outcome of this work: using the quark matter EoS derived from a covariant nonlocal 
chiral quark model, generalized here by including a vector interaction, we obtain hybrid star configurations which 
fulfill the modern constraints on high masses and radii of compact stars discussed in the Introduction. Moreover, the 
deconfinement transition for isospin symmetric matter occurs at about threefold nuclear saturation density, and results 
in a sufficient softening of the EoS thus circumventing a violation of the flow constraint from heavy-ion collisions. 
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FIG. 1: (Color online) Pressure as a function of the baryochemical potential for isospin symmetric matter. The nuclear matter 
phase is modelled by the DBHF equation of state with the Bonn- A potential [33j] (solid line) and the results for the covariant, 
nonlocal chiral quark model are given for different scaled vector coupling strengths g — 0.03, 0.05, 0.07, 0.09 (dotted, dash- 
dotted, dashed and dash-double-dotted curves, respectively) and scaled diquark coupling strengths of h = 0.70 (left panel) and 
h = 0.74 (right panel). A phase transition to quark matter is obtained at the crossing of nuclear and quark matter curves, for 
discussion see text. 
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FIG. 2: (Color online) Pressure as a function of the density for isospin symmetric matter. The phase transition to quark matter 
softens the EoS at densities above 0.55 fm -3 , relative to the nuclear DBHF EoS, thus allowing to fulfill the flow constraint 
derived in Ref. [13] . Line styles as in Fig. [T] 



IV. CONCLUSIONS 



We have generalized in this work a recently developed covariant, nonlocal chiral quark model by including a vector- 
vector four quark interaction which leads to a stiffening of the corresponding quark matter EoS. This allows us 
to describe compact stars possessing a quark matter core and being in agreement with the modern compact star 
phenomenology, which suggests stars with maximum masses of ~ 2 Mq and radii of 12 — 13 km. 
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FIG. 3: (Color online) Mass-radius relationships for neutron star configurations (DBHF EoS, solid line) and hybrid star 
configurations with hadronic shell (DBHF EoS) and color superconducting quark matter core (nonlocal chiral quark model 
EoS) for two parameter sets characterized by the coupling ratios h = H/Gs and g — Gv/Gs- Dashed line: h = 0.74, g = 0.07; 
dash-dotted line: h = 0.70, g = 0.05. 

We show that the traditional application of Gibbs conditions for phase equilibrium to construct a phase transition 
between deconfined quark matter - described in the new approach developed here - and nuclear matter - described 
by the DBHF approach - with the Bonn-A potential, is very sensitive to tiny changes of parameter values: quark 
and hadronic matter EsoS behave rather similarly in the vicinity of the phase transition and lead to a crossover-like 
behavior for the hybrid EoS. This reconfirms an earlier discussed 'masquerade' effect [1, HH, HU for hybrid stars within 
the present approach. After exploring the allowed range of model parameters, we find values for which both the flow 
constraint from heavy-ion collisions and the mass constraint for hybrid stars are satisfied. 

We want to point out that the deconfinement phase transition should be a rather robust phenomenon of the high- 
density EoS. Therefore, the mentioned 'masquerade' effect may point to a deficiency in the two phase approach for 
the description of hybrid EsoS. It is a demanding task to develop unified approaches for quark/nuclear matter on the 
basis of chiral quark models in which nucleons and mesons appear as relativistic bound states of quarks and antiquarks 
(for first steps in this direction, see Ref. [H, H(J). Under conditions of high density and/or temperature these bound 
states dissociate into continuum correlations (resonances) in quark matter within a Mott transition (see [5ll | for a 
model calculation within a nonrelativistic Green functions approach, which has also lead to an early prediction of 
stable quark matter cores in compact stars [Hj]). 

We believe that the covariant nonlocal chiral quark model presented here could be developed to a more elaborate 
approach unifying quark and hadronic matter descriptions on the quark level. The realization of such a project is 
beyond the scope of the present work. 
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APPENDIX A: DETAILS OF THE NONLOCAL MODEL FOR QUARK MATTER 

In this Appendix we show some explicit expressions corresponding to the nonlocal chiral quark model considered 
in Subsect. Ill Bl The determinant of the 48 x 48 matrix S' 1 appearing in Eq. ([7]) can be analytically calculated. In 
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this way one obtains 



QMFA = 



2G S 



A 2 



2H 2G V 



d 4 P 



4 E hai 2 , 



c = r,g,b 



where 



Ar = 



(Pic? + (Kc) 2 [Plf + (Kd 2 + (1 - 8k) A p2 A p 2 + 2p + ■ p- c + 2 W uc £* c 



with the following definitions: 



P% = if, P4 T ifJ-fc) , 
m f + a g{p^ c 2 ) , 



yP 



A 9 



(Al) 



(A2) 



(A3) 
(A4) 



where / = u,d, and c — r, g, b. Notice that due to the symmetry between red and green colors one has pf r = jif g . 

In general, for finite values of the current quark mass, Q MFA turns out to be divergent. We have used here a 
regularization procedure in which we add and subtract the thermodynamical potential for a free quark gas, namely 



r)MFA 

"(reg) 



QMFA 



Qfree 
( non r eg) 



ree 
eg)/c 



(A5) 



f,c 



In the right hand side, ^(™o nreg ) is obtained from Eq. (|Aip just by setting A = a — u> — 0, while for each fermion 
species i the regularized free thermodynamical potential ^fcg)t 1S §i ven by 



p>ficc 



24tt 2 



m\ F{ni/mi) , 



(A6) 



with 



F(x) = 2x{x 2 - 5/2) \Jx 2 - 1 + 3 In (x + \J x 2 - 1 ) 



(A7) 



The sum in Eq. (|A5|) extends over all quark flavors and colors. In the case of leptons, since they can be treated as 
free particles, the corresponding contribution to the compact star thermodynamical potential is simply given by 



free 
(reg)Z 



(A8) 



with I = e,fj,. 

Finally, according to our regularization prescription, the fermion densities appearing in Eq. () 1 3(1 are given by 



Pfc 



dfifc 



aofroc 
""(nonrog 



free 
H>eg)/ C ' 



_ „free 

Pi - P(reg); > 



(A9) 



where / = u,d, c = r, g, b, and I — e, fi. The fermion density p^eg) °^ a ^ ree P ar ticle gas can be easily obtained from 
the regularized thermodynamical potential in Eq. (|A6[) . yielding 



free 
P(rcg), 



1 

3~tt 



— 2 ip 2 ~ ™\? 12 



(A10) 



while explicit expressions for the partial derivatives in Eq. (|A9[) can be obtained from the results quoted in Ref. [32| . 
taking the T = limit. 
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